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-^We  present  a  brief  survey  of  the  current  state  of  the  art  in  network 
reliability.  We  survey  only  exact  methods;  Monte  Carlo  methods  are 
not  surveyed.  ^ 

^Most  network  reliability  problems  are,  in  the  worst  case,  NP-hard  and 
are,  in  a  sense,  more  difficult  than  many  standard  combinatorial  optimi¬ 
zation  problems.' 

.^Although  the  above  sounds  very  discouraging,  there  are  in  fact  linear 
and  polynomial  time  algorithms  for  network  reliability  problems  of 
special  structure. 

^We  review  general  methods  for  network  reliability  computation  and  dis¬ 
cuss  the  central  role  played  by  domination  theory  in  network  reliability 
computational  complexity.  We  also  point  out  the  connection  with  the 
more  general  problem  of  computing  the  reliability  of  coherent  structures  , 
(c.f.  IwdoiHmd  Pruschau  The  class  of  coherent  structures 

contains  both  directed  and  undirected  networks  as  well  as  logic  (or 
fault)  trees  without  not  gates.  This  is  a  rich  area  for  further  research. 
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A  SURVEY  OF  NETWORK  RELIABILITY 


by 

Avinash  Agrawal  and  Richard  E.  Barlow 


1.  INTRODUCTION 

Figure  1  is  a  well-known  example  of  a  computer  communication  network. 

A  two- terminal  reliability  problem  would  be  that  of  computing  the  probability 
that  the  distinguished  node  labelled  UCSB  can  communicate  with  the  distin¬ 
guished  node  labelled  CMU  via  some  set  of  arcs  or  edges.  Edges  may  be  sub¬ 
ject  to  failure.  In  this  paper,  the  edge  success  probabilities  are  assumed 
known  and  the  associated  success  events  are  assumed  independent  given  those 
probabilities.  A  typical  network  reliability  problem  is  to  calculate  effi¬ 
ciently  the  probability  that  a  specified  set  of  nodes  can  communicate  with 
each  other  at  a  given  time.  Edges  may  be  either  directed  or  undirected,  the 
direction,  when  it  applies,  will  be  denoted  by  an  arrow  on  the  edge. 

In  this  paper,  we  present  a  brief  survey  of  the  current  state  of  the 
art  in  network  reliability.  We  survey  only  exact  methods;  Monte  Carlo  methods 
are  not  surveyed.  There  are  many  papers  in  this  field,  some  give  methods  to 
simplify  or  solve  the  problems  while  others  calculate  the  complexity  of  net¬ 
work  reliability  problems.  In  1975,  Rosenthal  showed  that  certain  fault  tree 
and  network  reliability  problems  are  inherently  difficult  and  almost  certainly 
have  no  fast  algorithm.  Since  then,  this  class  has  grown  to  contain  a  number 
of  other  network  reliability  problems  and  some  of  the  questions  that  he  raised 
regarding  other  related  problems  have  been  answered. 


RAND  ALB  BBN  HARVARD  BURROUGHS 

ARPA  COMPUTER  NETWORK 


FIGURE  1 


We  begin  with  a  very  brief  introduction  to  the  theory  of  computational 
complexity  and  NP-completeness.  More  details  can  be  found  in  Garey  and 
Johnson  (1979).  An  algorithm  is  a  polynomial  time  algorithm  if  for  a 
problem  of  size  n  ,  its  running  time  is  bounded  by  a  polynomial  in  n  . 

Any  algorithm  that  is  not  a  polynomial  time  algorithm  is  generally  referred 
to  as  an  exponential  time  algorithm.  In  combinatorics,  the  so-called 
"satisfiability"  problem  is  in  the  class  NP-complete,  i.e.,  given  an  arbi¬ 
trary  Boolean  expression  in  product  of  sums  form,  determine  whether  or  not 
there  exists  an  assignment  of  values  TRUE  or  FALSE  to  the  variables  which 
makes  the  entire  expression  TRUE  (a  Boolean  expression  can  be  obtained  for 
every  network  in  terms  of  Boolean  indicators  for  the  edges  and  nodes) .  A 
problem  P  is  said  to  belong  to  the  class  NP-complete  if  (i)  given  a  pur¬ 
ported  solution  its  validity  can  be  checked  in  polynomial  time,  (ii)  the 
existence  of  an  algorithm  to  solve  P  in  polynomial  time  implies  the  exis¬ 
tence  of  algorithms  to  solve  the  satisfiability  problem  in  polynomial  time. 
It  is  generally  believed  that  no  polynomial  time  algorithm  exists  for  any 
of  the  NP-complete  problems.  Any  problem  not  NP-complete  but  which  can 
be  proved  to  be  at  least  as  hard  as  NP-complete  problems  is  known  as  an 
NP-hard  problem. 

Most  network  reliability  problems  are,  in  the  worst  case,  NP-hard 
[Ball  (1977  and  1980)].  Network  reliability  problems  are,  in  a  sense,  more 
difficult  than  many  standard  combinatorial  optimization  problems.  That  is, 
given  a  tentative  solution  to  a  combinatorial  problem,  often  its  correct¬ 
ness  can  be  determined  in  polynomial  time.  However,  given  a  purported  so¬ 
lution  to  a  reliability  problem,  it  cannot  even  be  checked  without  comput¬ 
ing  the  reliability  of  the  network  from  the  beginning. 
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Although  the  above  remarks  sound  very  discouraging,  there  are  in  fact 
linear  and  polynomial  time  algorithms  for  network  reliability  problems  of 
special  structure.  These  will  be  reviewed  in  Section  3.  In  Section  2, 
we  review  general  methods  for  network  reliability  computation  and  discuss 
the  central  role  played  by  domination  theory  in  network  reliability  com¬ 
putational  complexity.  We  also  point  out  the  connection  with  the  more  gen¬ 
eral  problem  of  computing  the  reliability  of  coherent  structures  [cf.  Barlow 
and  Proschan  (1981)].  The  class  of  coherent  structures  contains  both  di¬ 
rected  and  undirected  networks  as  well  as  logic  (or  fault)  trees  without 
not  gates  [cf.  R.  R.  Willie  (1978)].  This  is  a  rich  area  for  further  re¬ 
search. 


2.  COMPUTATIONAL  COMPLEXITY  OF  FACTORING  ALGORITHMS:  DOMINATION  THEORY 


Consider  a  graph  G  ■  (V,E)  with  vertex  set  V  and  edge  set  E  = 

{1,2,  n}  .  Vertices  do  not  fail,  but  at  an  instant  of  interest,  an 

edge  i  has  reliability  p^  ,  independent  of  the  states  of  other  edges. 
Figure  2  is  a  network  graph  with  edge  set  E  -  {l,2,3,4,e}  .  Since  no  ar¬ 
rows  appear  on  edges,  this  is  an  undirected  network.  The  graph,  together 
with  the  distinguished  nodes  K  ■  {s,t}  ,  define  a  network  reliability 
problem  -  namely,  calculate  the  probability  that  s  and  t  can  communi¬ 
cate.  The  distinguished  nodes  also  define  the  topology  of  the  problem  - 
namely,  the  family  of  minimal  path  sets,  P  «  [{l,3},{l,2,4},{e,4},{e,2,3}]  . 
(A  minimal  path  set  is  a  minimal  set  of  elements  whose  functioning  implies 
that  the  distinguished  nodes  can  communicate.)  A  different  set  of  distin¬ 
guished  nodes  would  define  a  different  topology. 


UNDIRECTED  TWO-TERMINAL  NETWORK 
FIGURE  2 


Figure  3  is  a  rooted  directed  graph  with  distinguished  node  set  K  ■ 
{s,t,v}  .  The  "root”  is  vertex  s  ,  the  unique  node  with  no  entering  edges. 
The  problem  is  to  compute  the  probability  that  s  can  communicate  with 
both  v  and  t  .  The  topology  for  this  problem  is  defined  by  the  family 
of  minimal  path  sets  P  -  [{2, 5), {1,3, 5}, {2, e, 4}, {1,3, 4}, {1,2, 4}]  . 
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DIRECTED  NETWORK 
FIGURE  3 

Graphically,  a  minimal  path  for  this  problem  is  a  "rooted  tree."  For 
example,  the  following  acyclic  rooted  directed  graph  is  a  "tree"  correspond¬ 
ing  to  Figure  3. 


v 


ROOTED  TREE 

FIGURE  4  ' 

Since,  historically,  not  much  emphasis  has  been  placed  in  the  reli¬ 
ability  literature  on  computational  complexity  and  because  these  problems 
are,  in  general,  NP-hard,  many  different  algorithms  have  been  suggested  to 
solve  these  problems.  Many  algorithms  are  based  on  minimal  path  sets  and/or 
cut  sets  (a  minimal  cut  set  is  a  minimal  set  of  elements  whose  failure  im¬ 
plies  that  some  distinguished  nodes  cannot  communicate) .  In  general ,  how¬ 
ever,  it  is  neither  necessary  nor  desirable  to  find  the  family  of  minimal 
paths  or  cut  sets  in  order  to  calculate  network  reliability. 
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The  reliability  of  a  graph,  G  ,  with  distinguished  nodes  K  ,  is  the 
probability,  R^G)  ,  that  all  elements  of  at  least  one  minimal  path  set 
are  working  or  one  minus  the  probability  that  all  elements  of  at  least  one 
minimal  cut  set  have  failed.  Note  that  R^G)  depends  on  the  distinguished 
node  set  K  C  V  as  well  as  G  .  Different  methods  exist  to  evaluate  this 
probability.  These  methods  are  quite  general  and  can  be  used  relative  to 
any  system  reliability  problem.  We  shall  roughly  classify  them  into  three 
main  classes.  Let  denote  the  event  that  all  elements  in  the  i^  mini¬ 
mal  path  set  are  functional  and  denote  the  complement  of  this  event. 

Let  p  be  the  number  of  minimal  path  sets. 


(i)  The  Inclusion-Exclusion  Method 


Rjf (G)  -  P  U  A.  -  l  P(A  )  -  l  l  P[A  A  ] 
|_i-l  i=l  i=l  j<i  j 


+  . . . 


+  <-ir  ...  Ap]  . 


(2.1) 


If  there  are  p  path  sets,  then  this  calculation  involves 
2P  -  1  terms.  In  some  cases,  two  different  intersections  of 
A^'s  will  correspond  to  the  same  event  so  that  these  different 
intersections  of  A^'s  will  have  the  same  probability.  If  one 
Intersection  consists  of  an  odd  number  of  A^'s  and  another 
intersection  consists  of  an  even  number  of  A^'s  ,  their  prob¬ 
abilities  in  the  inclusion-exclusion  expression  will  cancel. 
Satyanarayana  and  Prabhakar  (1978)  and  Satyanarayana  (1982) 
give  algorithms  that  generate  only  the  non-cancelling  terms. 

In  the  reduced  inclusion-exclusion  expression,  there  will  be  a 
coefficient  for  the  term  corresponding  to  the  event  AjA^  . . .  A^ 
(This  event  will  in  general  correspond  to  several  different 
intersections.)  The  coefficient  for  this  term  is  called  the 


* .  wRu  .9.  .1 ..  sRw. .  1 1  .,.^^1 . .  .9 .  j—.  ... 
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signed  domination,  d^(G)  .  This  coefficient  is  the  number  of 
odd  formations  of  G  minus  the  number  of  even  formations  of 
G  .  (A  formation  is  a  set  of  minimal  path  sets  whose  union 
constitutes  the  edge  set  E  .  It  is  odd  (even)  if  it  has  an  odd 
(even)  number  of  minimal  path  sets  in  it.)  The  absolute  value  of 
this  coefficient  is  called  the  domination,  D^(G)  =  | dR (G) |  .  As 
we  shall  see,  this  number  is  a  measure  of  the  computational  com¬ 
plexity  of  certain  factoring  algorithms  for  undirected  networks. 

The  network  corresponding  to  Figure  2  has  6  formations  of  which 
4  are  odd  and  2  are  even  so  that  d^(G)  =  2  in  this  case. 

(ii)  Sum  of  Disjoint  Products 

^(G)  -  P(A1)  +  P(IlA2)  +  ...  +  P(AjA2  ...  Ap_jAp)  .  (2.2) 

For  p  path  sets,  there  are  p  terms  but  the  time  needed  to 
generate  each  term  may  be  exponential  in  p  .  Most  methods  based 
on  Boolean  techniques  belong  to  this  class.  Such  methods  have 
been  proposed  by  Abraham  (1979) ,  Fratta  and  Montanari  (1973) , 

Aggarwal,  Mlsra  and  Gupta  (1975),  etc.  One  proposed  by  Lee  (1979) 

for  flow  networks  can  easily  be  modified  by  using  concepts  of 

Satyanarayana  (1982)  to  give  the  reliability  of  communication 

networks.  This  method,  based  on  the  idea  of  using  backtracking 

and  depth  first  search  introduced  by  Gabow  and  Myers  (1978) ,  is  j 

the  most  efficient  method  in  this  class  for  the  undirected  all-  j 

terminal  and  the  directed  source-to-all-terminal  problems;  i.e.,  i 

i 

the  distinguished  set  of  vertices  is  the  set  of  all  vertices  of  the  { 

graph.  The  amount  of  computational  work  required  is  proportional  to  1 

the  number  of  spanning  trees  or  rooted  spanning  arboresences.  [See  i 

I 

Deo  (1974)  for  a  detailed  discussion  of  graph  theory  concepts.]  \ 

\ 


.  *  J  '  -L' 
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(iii)  Pivotal  Decomposition  or  Factoring 

If  R^(G  I  ®)  is  the  reliability  of  G  under  the  condition 
that  edge  e  is  working  and  R^G  |  e)  the  reliability  of  G 
under  the  condition  that  edge  e  is  not  working,  then  using 
the  pivotal  decomposition  [cf.  Barlow  and  Proschan  (1981)], 

R^G)  -  p^CG  |  e)  +  (1  -  p^R^G  j  i)  .  (2.3) 

R^G)  for  any  graph  G  can  be  computed  by  repeated  applica¬ 
tion  of  this  decomposition.  Undirected  graphs  have  some  special 
properties  that  can  be  used  to  simplify  this  method.  If  the 
vertices  are  assumed  to  be  working,  then  R^(G  |  e)  is  the  same 
as  Rj j(Gg)  where  Gg  is  the  graph  obtained  from  G  by  delet¬ 
ing  edge  e  and  merging  its  end  points.  Similarly,  R^G  |  e) 
is  the  same  as  RR(G  -  e)  where  G  -  e  is  the  graph  with  e 
deleted  but  no  vertex  deleted.  These  schemes  have  been  dis¬ 
cussed  by  a  number  of  writers,  among  them,  Hoskowitz  (1958)  and 
Misra  (1970).  In  the  literature,  the  above  relationship  has 
been  referred  to  as  the  Factoring  Theorem .  It  is  important  to 
note  that  this  method  can  be  employed  using  the  graph  repre¬ 
sentation  but  without  knowing  the  minimal  path  sets.  However, 
unless  some  kind  of  probability  reductions  are  performed  (e.g., 
parallel  and  series  reductions)  after  each  pivot,  this  method 
will  be  equivalent  to  state  space  enumeration. 

The  Factoring  Algorithm  for  Undirected  Graphs 
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reductions  as  the  algorithm  proceeds.  We  recursively  apply  Equation  (2.3), 
making  probability  reductions  as  we  go.  The  edge  selection  strategy  for 
computing  (2.3)  attempts  to  avoid  creating  subgraphs  with  irrelevant  edges. 

Figure  5  shows  the  binary  computational  tree  resulting  from  using  the 
factoring  algorithm  for  the  two-terminal  problem  corresponding  to  the  graph 
at  the  top  of  Figure  5  with  distinguished  nodes  s  and  t  ,  so  that  K  * 
{s,t)  .  In  the  initial  step,  we  pivot  on  edge  e  forming  two  subgraphs: 

Gg  ,  corresponding  to  e  working,  and  G  -  e,  corresponding  to  e  failed. 
Parallel  and  series  probability  reductions  are  now  possible.  In  a  parallel 
probability  reduction,  two  edges,  say  2  and  3  in  Gg  ,  are  replaced  by  a 
single  edge  with  associated  probability  P2  +  P3  ~  P2P3  *  Likewise  in  Gg  , 
this  new  edge  with  edge  5  are  in  series.  The  two  edges  in  series  are  then 

replaced  by  a  single  edge  with  associated  reliability  (p£  +  p^  -  P2p3)?5  • 

(In  Figure  5,  the  parallel  and  series  reductions  are  not  shown.)  Pivoting 
now  proceeds  on  edge  4  resulting  in  two  additional  subgraphs,  each  of  which 
can  be  reduced  to  a  single  edge  by  series  and  parallel  probability  reduc¬ 
tions.  The  "leaves"  of  the  tree  are  the  four  subgraphs  at  the  bottom  of 
the  tree. 

If  each  edge  i  has  probability  p  of  working,  it  can  be  shown 
using  the  binary  computational  tree  (illustrated  in  Figure  5)  that  the  sys¬ 
tem  reliability  in  this  case  is 

R^G)  -  p2(((((p  it  p)p)  a  p)p)  11  p)  +  p(l  -  p)(((p  a  p)p)(p2  a  p)) 

+  p(l  -  p)((p(p  a  p))  a  (p2))p  +  (1  -  p)2((p3  a  p)p2) 

■  (p3  +  +  p3  -  5p^  +  4p^  -  p®)  +  (2p^  -  p3  -  4p^  +  4p2  -  p®) 

4  .5  6,-7  8..  ,,  3  -4,-5  ,6.-7  8. 

+  (3p  -  4p  -  p  +  3p  -  p  )  +  (p  -  2p  +  2p  -  3p  +  3p  -  p  ) 

-  3  .  ,  4  ,5  i/7  .8 

■  2p  +  4p  -  2p  -  13p  +  I4p  -  4p  . 


The  lower  case  "ip"  operator,  u  ,  corresponds  to  calculating  the  reli¬ 
ability  of  parallel  edges;  i.e., 

"i  *  ■  »i  +  -  Vj  • 

In  Figure  5,  p^  ■  1  -  p^  .  The  four  graphs  at  the  bottom  of  the  tree 
are  the  leaves  of  the  tree  and  each  has  domination  one  since  each  is 
series-parallel  reducible  (Chang,  1981).  The  domination  of  the  top  graph, 
it  turns  out,  is  D^(G)  *  4  (the  number  of  "leaves"  at  the  bottom  of  the 
tree)  and  the  tree  has  2D^(G)  -  1  ■  7  nodes  so  that  the  computational  run¬ 
ning  time  is  proportional  to  the  domination.  Satyanarayana  and  Chang  (1983) 
found  that  in  general  the  number  of  leaves  in  the  binary  computational  tree 
using  a  factoring  algorithm  with  series  and  parallel  probability  reductions 
is  at  least  equal  to  the  domination.  Using  a  simple  edge  selection  strate¬ 
gy,  they  further  showed  that  it  is  possible  to  create  a  backtrack  structure 
which  has  exactly  DR(G)  leaves  where  DR(G)  is  the  domination  of  G  . 
Therefore,  this  edge  selection  strategy  is  optimal  for  factoring  algorithms 
using  series  and  parallel  probability  reductions. 

Domination  Theory  for  Coherent  Systems 

A  network  graph  G  with  distinguished  nodes  K  C  V  has  a  topology 
defined  by  the  minimal  path  sets  P  *  [P^ »I>2 *  •••»  Ppl  •  ^*e  network  may 

have  both  directed  and  undirected  edges.  By  definition,  all  nodes  K  can 
communicate  if  and  only  if  all  edges  in  at  least  one  minimal  path  set  oper¬ 
ate.  In  Barlow  and  Proschan  (1981),  a  set  of  edges  or  components  E  ,  and 
a  family  of  minimal  path  sets,  P  ,  is  called  a  coherent  system,  (E,P)  , 

if  P  ■  [P^^,  ...,  P  ]  is  a  minimal  family,  i.e.,  no  P^  is  contained 

,  P 

in  another  member  of  the  family  and  E  ■  UP. 

i-1  1 


Coherent  systems  include 
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all  network  graphs  as  well  as  logic  trees  (or  fault  trees)  without  not  gates. 
A  k-out-of-n  system  with  n  >  2  and  1  <  k  <  n  Is  a  coherent  system  which 
cannot  be  represented  as  a  network  (unless  replicated  edges  are  allowed) . 

The  two-terminal  network  In  Figure  2  will  be  used  to  illustrate  Ideas. 
For  this  example,  E  -  {l,2,3,4,e}  while  P  *  [{1, 3}, {1,2, 4), {e, 4}, (e, 2,3}] 
By  pivoting  on  component  e  e  E  ,  we  create  two  subsystems,  correspond- 


v". 

’*/-■ 

ing  to  the  system  with  e  failed  and  to  the 

system  with  e  perfect,  res- 

1 

pectively.  Let 

P(e)  -  [P±  |  a  e  P±  and  ?± 

e  P) 

and  P(e*)  -  [P±  |  e  t  P± 

?: 

and  Pi  e  P]  . 

Then 

r.- 

r<r 

P  -  P(e)  U  P(e')  . 

I 

In  our  example. 

P(e)  -  [{e,4},{e,2,3}]  and 

P(e') 

-  [{1,3}, {1,2, 4}]  .  In 

all  cases,  U  P^  C  E  -  e  .  In  our  example,  U 

P^Pte')  P±eP(e') 


E  -  e  so  that 


(E  -  e,P(e'))  is  coherent  and  corresponds  to  our  system  with  e  failed. 

If  UPi  C  E-e  ,  then  (E  -  e,P(e'))  would  have  no  formations  so  that 
P ±eP(e*) 

in  this  case  d(E  -  e,P(e'))  *  0  and  (E  -  e,P(e'))  would  be  noncoherent. 
(In  this  setup,  d(E,P)  denotes  signed  domination.) 

To  describe  a  system  with  e  perfect,  let 


% 

& 


5: 


P  —  e  ■  [P^  —  e,?2  -  e ,  •••»  Pp  ~  ®]  • 


If  e  i  P^  ,  then  P^  is  included  as  it  is.  Let  M[P  -  e]  be  the  set 
minimization  of  P  -  e  .  In  our  example, 

P  -  e  -  [{1,3}, {1,2, 4}, {4}, {2, 3}] 

and 


M[P  -  e]  -  [{1,3}, {4}, {2, 3}] 


since  {4}  C  {1,2,4}  .  In  this  case. 


U  A±  -  E  -  e  , 

AieM[P-e] 

so  that  (E  -  e,M[P  -  e])  corresponds  to  our  example  system  with  e 
perfect.  In  general,  we  only  know  that 

U  A±  CE-e 
AiEM[P-e] 

so  that  (E  -  e,M[P  -  e])  might  be  noncoherent. 

The  following  signed  domination  theorem  is  proved  in  Barlow  (1982) . 

It  was  first  proved  by  Satyanarayana  and  Chang  (1983)  for  undirected  networks. 

Theorem  2.0:  (Signed  Domination  Theorem) 

For  any  coherent  system  (E,P)  and  e  e  E  , 

d(E,P)  -  d(E  -  e,M[P  -  e])  -  d(E  -  e,P(e’))  . 

In  our  example,  d(E  -  e,M(P  -  ej)  -  1  and  d(E  -  e,P(e'))  -  -1 
so  that  d(E,P)  -  2  . 

Using  Theorem  2.0  and  Induction,  it  is  easy  to  verify  the  following 
corollary  for  undirected  networks. 

Corollary: 

For  a  coherent  system  corresponding  to  a  K-terminal  undirected  network 
problem, 

d(E,P)  -  (-l)n’v+1D(E,P) 

where  D(E,P)  ■  | d (E , P) |  is  the  domination,  n  is  the  number  of  edges 
and  v  is  the  number  of  vertices  or  nodes. 
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*'  Satyanarayana  and  Chang  (1983)  proved  the  following  domination  theorem. 

t,  Theorem  2.1:  (Domination  Theorem) 

For  any  K-terminal  undirected  graph, 

D(E,P)  -  D(E  -  e,M[P  -  e])  +  D(E  -  e,P(e’))  . 

Proof : 

By  Theorem  2.0, 

d(E,P)  -  d(E  -  e,M[P  -  e])  -  d(E  -  e,P(e’))  . 

By  the  corollary  to  Theorem  2.0, 

d(E,P)  -  (-l)n"v+1D(E,P) 

i 

if  the  undirected  graph  has  n  edges  and  v  nodes.  The  subgraph  corre- 
.<■  sponding  to  e  working  has  signed  domination 

d(E  -  e,M[P  -  e])  -  (-1) (n"1)"(v"1)+1D(E  -  e,M[P  -  e)) 

since  e  has  been  contracted  and  its  two  end  vertices  merged  into  one. 
Similarly,  the  subgraph  corresponding  to  e  failed  has  signed  domination 

d(E  -  e.P(e’))  -  (-1) (n-1)”v+1D(E  -  e,P(e’)) 

since  only  edge  e  is  deleted.  Hence, by  Theorem  2.0, 

d(E,P)  -  (-l)(n"1)"(v“1)+1D(E  -  e,M[P  -  e]) 

-  (-l)(n_1)"V+1D(E  -  e,P(e')) 

*  -  (-l)(n“1)_(v"1)+1[D(E  -  e,M[P  -  e])  +  D(E  -  e,P(e’))]  .  Q.E.D. 


The  domination  theorem  is  in  general  not  true  for  directed  graphs. 

For  example,  the  domination  of  the  cyclic  directed  graph  in  Figure  3  is  zero. 
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(In  fact,  all  directed  networks  with  cycles  have  domination  zero  [cf.  R.  R. 
Willie  (1980)].)  However,  if  we  pivot  on  e  ,  we  obtain  a  subgraph  corre¬ 
sponding  to  e  deleted  which  has  positive  domination.  Hence,  the  domina¬ 
tion  theorem  is  not  true  for  this  example. 

From  Theorem  2.1,  it  can  be  shown  that  the  computational  complexity  of 
any  factoring  algorithm  for  undirected  networks  based  on  pivoting  and  par¬ 
allel  and  series  reductions  will  require  at  least  2D(E,P)  -  1  steps. 

(Every  binary  computational  tree  has  number  of  nodes  equal  to  twice  the 
number  of  leaves  minus  one.)  By  using  an  edge  selection  strategy  for  which 
both  D(E  -  e,M[P  -  e])  >  0  and  D(E  -  e,P(e’))  >  0  ,  the  number  of  steps 
is  exactly  2D(E,P)  -  1  .  This  is  illustrated  by  Figure  5. 

The  factoring  algorithm ,  namely  pivoting  followed  by  parallel  and 
series  probability  reductions  and  repeated  until  only  single  edges  or  "K 
trees"  are  obtained,  can  be  applied  to  any  coherent  system  reliability 
problem.  (A  K-tree  of  a  graph  G  with  respect  to  K  is  any  minimal  graph 
which  connects  all  the  distinguished  nodes  in  K  .)  However,  the  optimal 
edge  selection  strategy  is  only  known  because  of  the  domination  theorem. 

The  domination  itself  is  mainly  of  theoretical  interest.  In  the  worst  case, 
the  domination  will  be  exponential  in  the  number  of  edges  and  would  normally 
not  be  computed.  It  does,  however,  offer  a  sort  of  theoretical  benchmark. 
For  example,  if  we  delete  all  arrows  from  a  directed  graph,  the  resulting 
undirected  graph  will  have  a  domination  value  which  will  provide  an  upper 
bound  on  the  best  factoring  algorithm  for  the  directed  network.  This  is  so 
because  the  domination  of  the  two  subgraphs  created  by  pivoting  on  an  edge 
in  the  undirected  graph  will  be  at  least  as  much  as  the  domination  of  the 
undirected  graphs  corresponding  to  the  subgraphs  obtained  by  pivoting  on 
the  same  edge  in  the  directed  graph. 
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3.  SPECIAL  STRUCTURES 

Although  the  factoring  algorithm  can.  In  principle,  solve  all  reli¬ 
ability  problems,  It  is.  In  the  worst  case,  an  exponential  time  algorithm. 
For  very  large  networks,  we  need  linear  or  polynomial  time  algorithms  in 
order  to  calculate  system  reliability  in  "reasonable"  computing  time.  By 
introducing  additional  probability  reductions,  such  algorithms  have  been 
found  for  both  directed  and  undirected  network  graphs  of  special  structure. 

An  undirected  graph  G  ■  (V,E)  is  said  to  be  basically  series-parallel 
if  the  graph  (without  distinguished  nodes)  can  be  reduced  to  a  single  edge 
by  series  and  parallel  replacements.  A  replacement  as  opposed  to  a  prob¬ 
ability  reduction  does  not  involve  the  probability  measure  which  may  be 
associaced  with  the  graph.  For  example,  Figure  2  can  be  reduced  to  a  single 
edge  by  series  and  parallel  replacements  as  follows: 


SERIES  AND  PARALLEL  REPLACEMENTS 
FIGURE  6 

Note  that  vertices  s  and  t  are  no  longer  distinguished.  Edges  e  and 
1  are  replaced  by  a  dotted  line  using  a  series  replacement  as  are  edges  3 
and  4.  Finally,  the  remaining  three  edges  in  parallel  are  replaced  by  a 
single  edge.  No  probability  calculations  are  involved. 


The  network  in  Figure  1  on  the  other  hand  is  not  basically  series- 
yarallel.  Neither  is  the  network  in  Figure  7. 


EXAMPLE  NETWORK  WHICH  IS  NOT  BASICALLY  SERIES-PARALLEL 

FIGURE  7 

A  directed  network  Is  basically  series-parallel  if  the  underlying  un¬ 
directed  graph  (without  arrows  on  edges)  is  basically  series-parallel. 

For  example,  the  network  in  Figure  3  is  basically  series-parallel. 

A.  Satyanarayana  and  R.  K.  Wood  (1982)  provide  linear  time  algorithms 
for  calculating  the  K  terminal  reliability  of  undirected  networks  which  are 
basically  series-parallel.  They  introduce  probability  reductions  called 
polygon-to-chain  reductions  to  accomplish  this. 

A.  Agrawal  and  A.  Satyanarayana  (1983)  provide  linear  time  algorithms 
for  calculating  the  source  to  K  terminal  reliability  of  rootedt  directed 
networks  which  are  basically  series-parallel.  One  node  in  K  is  designated 
the  root  and  the  reliability  problem  is  to  calculate  the  probability  that  the 
root  can  conmmicate  with  the  remaining  K  C  V  vertices. 
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The  Minimal  Domination  of  Undirected  Graphs 

The  minimal  domination  M(G)  of  a  graph  G  is  defined  by 

M(G)  “  minimum  DV(G) 

K:|K|-2 

where  K  is  a  distinguished  set  of  nodes  of  G  .  While  the  domination 
Dg(G)  depends  on  both  the  graph  G  and  the  distinguished  set  of  nodes 
K  ,  M(G)  obviously  depends  only  on  G  .  Whereas  D^(G)  =1  if  and  only 
if  G  is  reducible  to  a  "K-tree"  by  series  and  parallel  probability  re¬ 
ductions,  [Satyanarayana  and  Chang  (1983)],  M(G)  >1  if  and  only  if  G 
is  basically  series-parallel  [R.  K.  Wood  (1982)].  Thus,  the  graph  at  the 
top  of  Figure  5  has  DR(G)  ■  4  where  K  -  (s,t)  ,  but  M(G)  ■  1  since  G 
is  basically  series-parallel.  Using  series  and  parallel  probability  reduc¬ 
tions  and  the  polygon-to-chain  reduction  in  Table  1,  ^(G)  can,  in  this 

case,  be  computed  without  pivoting  so  that  a  linear  time  algorithm  exists 
for  this  problem  and  in  fact  for  all  such  problems  where  G  is  basically 
series-parallel .  There  are  6  additional  polygon-to-chain  reductions  nec¬ 
essary  to  treat  cases  where  |k|  >  2  .  [See  R.  K.  Wood  (1982).] 

Figure  8  provides  an  example  graph  where  the  domination  is  D^(G)  * 


FIGURE  8 
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''  The  minimum  domination  theorem  was  proved  by  R.  Procesi-Ciampi  (1981) . 

[See  also  Satyanarayana  and  Procesi-Ciampi  (1981).] 

% 

Theorem  3.1:  (Minimum  Domination  Theorem) 

For  any  undirected  graph  G  ■  (V,E)  , 

M(G)  -  M(Ge)  +  M(G  -  e)  . 

R.  K.  Wood  (1982)  used  this  and  other  properties  of  minimum  domination  to 
evaluate  the  computational  complexity  of  undirected  networks  relative  to 
pivoting  and  polygon-to-chain  reductions. 
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4.  CONCLUSION 

A  theoretical  breakthrough  occurred  in  1978  with  the  publication  by 
Satyanarayana  and  Prabhakar  which  first  introduced  the  idea  of  domination 
into  network  reliability.  This  paper  was  concerned  with  directed  graphs 
and  suggested  that  directed  cyclic  graphs  have  domination  zero,  a  result 
rigorously  proved  by  R.  R.  Willie  (1980).  In  1981,  Satyanarayana  and 
Chang  (see  their  1983  paper)  first  noticed  the  connection  between  domina¬ 
tion  theory  and  the  computational  complexity  of  factoring  algorithms.  This 
was  followed  by  the  notion  of  minimal  domination  and  the  development  of 
linear  time  algorithms  for  special  structure  graphs  [cf.  R.  K.  Wood  (1982)]. 

The  probability  measure  associated  with  the  graphs  in  the  literature 
surveyed  is  extremely  simplistic  and  unrealistic  for  many  practical  network 
reliability  problems.  However,  a  theoretical  basis  now  exists  for  more 
research  into  realistic  probability  measures.  At  the  present  time,  there  is 
no  comparable  theoretical  basis  for  analyzing  the  computational  complexity 
of  logic  (or  fault)  trees.  These  structures  are  perhaps  more  useful  in 
practice  than  networks.  Further  research  is  badly  needed  in  this  field. 
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